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d , Abstract. We consider the possibility to observationally differentiate the Standard Model 

I (SM) Higgs driven inflation with non-minimal coupling to gravity from other variants of 

CN ' SM Higgs inflation based on the scalar field theories with non-canonical kinetic term such 

, J^ . as Galileon-like kinetic term and kinetic term with non-minimal derivative coupling to the 

^T) ', Einstein tensor. 

■^ ' In order to ensure consistent results, we study the SM Higgs inflation variants by using the 

; . same method, computing the full dynamics of the background and perturbations of the Higgs 

tlT' I field during inflation at quantum level. 

,_^ ' Assuming that all the SM Higgs inflation variants are consistent theories, we use the MCMC 

^-H . technique to derive constraints on the inflationnoary parameters and the Higgs boson mass 

^ ! from their fit to WMAP7-|-SN-|-BA0 data set. We conclude that a combination of a Higgs 

l> I mass measurement by the LHC and accurate determination by the PLANCK satellite of the 

j_j ■ spectral index of curvature perturbations and tensor-to-scalar ratio will enable to distinguish 

among these models. We also show that the consistency relations of the SM Higgs inflation 

variants are distinct enough to differentiate the models. 
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1 Introduction 

There has been much interest recently in models of inflation in which the Standard Model 
(SM) Higgs boson non-minimally coupled to the Ricci scalar can give rise to inflation without 
the need for additional degrees of freedom to the SM [1-5]. This scenario is based on the ob- 
servation that the problem of the very small value of the Higgs quadratic coupling required 
by the Cosmic Microwave Background (CMB) anisotropy data can be solved if the Higgs 
inflaton has a large non-minimally coupling with gravity [6]. The resultant Higgs inflaton 
effective potential in the inflationary domain is effectively flat and can result in successful 
inflation for values of non-minimally coupling constant C, ~ 10^ (^-inflation), allowing for 
cosmological values for Higgs boson mass in a window in which the electroweak vacuum is 
stable and therefore sensitive to the field fluctuations during the early stages of the universe 

[7]. _ 

Limits of validity of (^-inflation have been recently debated by several authors. Specifically, 
Barbon & Espinosa (2009) argued that the large coupling of Higgs inflaton to the Ricci 
scalar makes this model invalid beyond the ultraviolet cutoff scale A^ ~ Mp/C, (hereafter 
Mp = 2.4 X 10^^ GeV is the reduced Planck mass) which is below the Higgs field expectation 
value at J\f e-foldings during inflation, (j) ~ ^/NMp/^/Q. As consequence, at the ultraviolet 
cutoff scale A^ at least one of the cross-sections of different scattering processes hits the uni- 
tarity bound [9] . The fact that the quantum corrections due to the strong coupling to gravity 
makes the perturbative analysis to break down at energy scales above A^ was interpreted 



as a signature of a new physics, implying higher dimensional operators at energies above 
A^. Few authors addressed the issue of ^-inflation naturalness with respect to perturbativity 
and unitarity violation in Jordan and Einstein frames [10-13]. They show that the apparent 
breakdown of the theory in the Jordan frame does not imply new physics, but a failure of 
the perturbation theory in the Jordan frame as a calculation method. These works demon- 
strate that, for inflation based on a single scalar field with large non-minimal coupling, the 
quantum corrections at high energy scales are small, making the perturbative analysis valid 
and consequently there is no a breakdown of unitarity at the energy scale A^. 
Recent works [14-16] revisit the issue of self-consistency of (^-inflation model emphasizing 
that the scale of unitarity violation depends on the size of the Higgs background field. They 
found the cutoff scale A^ higher than the relevant dynamical scales throughout the whole 
history of the universe, including the inationary epoch and reheating, indicating that the 
theory does not violate unitarity at any time. 

Although the above conclusion theoretically may allow to work within the regime of validity 
of the effective theory for inflationary calculations, other variants of Higgs inflation within the 
SM have been proposed. These models are based on scalar fleld theories with non-canonical 
kinetic term [17] and are motivated by some particle physics models of inflation such us the 
Dirac-Born-Infed inflation [18, 19]. 

A model with non-minimal derivative couplings was proposed in [20-22] in context of 
inationary cosmology, and recently, the inflation driven by a scalar field with a Higgs potential 
and non- minimal kinetic couplings to itself and to the curvature was proposed in [23, 24]. In 
the new model of Higgs inflation (E-inflation) an enhanced kinetic term in the Lagrangian 
propagating no more degrees of freedom than the minimally coupled scalar field in General 
Relativity (GR) is obtained by considering non-minimal derivative coupling of the SM Higgs 
boson to the Einstein tensor [25]. The power counting analysis indicates that the scale at 
which the unitarity bound is violated in this model is A(if) ~ {2H'^kY'^, where H is the 
Hubble expansion during inflation. Requiring the unitarity constraint i/ ^ A to be satisfied, 
it is shown that the non-minimal derivative coupling allows a Higgs boson self-coupling value A 
within the limits expected from the collider experiments while the cosmological perturbations 
in this model are consistent with present cosmological data [26, 27]. 

The Galileon models of inflation (G-inflation) are constructed by introducing a scalar 
fleld with a self-interaction whose Lagrangian is invariant under Galileon symmetry d^(j) ~^ 
dfj^(j) + bfj,, which maintains the equation of motion of the scalar fleld as a second-order differ- 
ential equation, preventing the theory from exhibiting new degrees of freedom and avoiding 
ghost or instability problems [28-30]. Recently, an inflation model in which inflation is 
driven by the SM Higgs scalar fleld with a Galileon- like kinetic term has been proposed [31]. 
The dynamics of background and perturbations of G-inflaton show new features brought by 
the modified kinetic term in the Lagrangian, when compared with the standard slow-roll 
inflation, such as: violation of the standard consistency relation [31, 32], violation of the 
Null Energy Condition [33], generation of isocurvature perturbations and large primordial 
non-Gaussianity [34-37]. 

Li this paper we analyze to possibility to distinguish among the SM Higgs inflation 
variants by using the existing cosmological and astrophysical measurements. We will adopt a 
similar philosophy as in [38], considering that the SM Higgs inflation variants are consistent 
theories and therefore should be subject to rigorous tests against experimental data. 
Cosmological constraints on the SM Higgs driven ("-inflation have been discussed by a number 



of authors [2-4]. These constraints are based on mapping between the RenormaUzation Group 
(RG) flow equations and the spectral index of the curvature perturbations obtained by the 
WMAP team, parametrized in terms of the number of e-foldings till the end of inflation. 
In a recent paper [39] we studied the implications of the full dynamics of background and 
perturbations of the Higgs field during ("-inflation at quantum level and derive constaints on 
the inflationary parameters and the Higgs boson mass from the analysis of WMAP 7-year 
CMB measurements complemented with astrophysical distance measurements [45, 46]. 
Our goal in this paper is to refine the computation for ^-inflation by considering in the analysis 
the improved RG flow equations, accounting at the same time for the possible uncertainties in 
the theoretical determination of the reheating temperature. We also make a similar analysis 
for G-inflation and E-inflation models and derive constraints on the inflationary parameters 
and Higgs boson mass from their fit to the same data set, ensuring in this way that the 
differences in the predictions of the observables are due to the differences in the underlying 
theories of the SM Higgs infiation and therefore can be used to distinguish among them. 

The paper is organized as follows. In Section 2 we introduce the there models of the 
SM Higgs inflation that we consider: (^-inflation [1], G-inflation [31] and E-inflation [25, 26]. 
In Section 3 we compute the RG improved Higgs field potential and in Section 4 we present 
our main results. In Section 5 we draw our conclusions. 

Throughout the paper we consider an homogeneous and isotropic flat background described 
by the Friedmann-Robertson- Walker (FRW) metric: 

ds^ = g^^^ (\x^ dx"" = -dt^ + o^(t)dx^ , (1.1) 

where a is the cosmological scale factor, k^ = SttMJ (where Mpi ~ 1.22 x lO"*^^ GeV is the 

Planck mass), overdots denotes the time derivatives and ^x = d/dX. 

We also remind that the canonical tree-level action describing the SM Higgs inflation is given 

by: 

Sq= j d^ Xy/^ 



— ^R -I- CsM 



(1.2) 



where R is the Ricci scalar and Csm is the tree-level SM Higgs Lagrangian: 

CsM = g^^v^n^v.n - A (n^n -v^f . (1.3) 

In the above expression Ti is the Higgs boson doublet, D^ is the covariant derivative with 
respect to SU{2) x SU{1) and v is the vacuum expectation value (vev) of the Higgs broken 
phase of the SM. Rotating the Higgs doublet such that Ti^ = (l/-v/2)(0, v + ({>) and assuming 
that (p is much greather than vev during inflation we have: 

where A is the tree-level Higgs self coupling and V{(j)) is the corresponding Higgs potential. 



CsM = Pi^,X)=X-Vicl)), x = --g^''d^d,cp, V{^) = -^\ (1.4) 



2 Models: background equations and cosmological perturbations 

In this section we derive the background equations of motion of the scalar fleld and present 
the power spectra of scalar and tensor perturbations derived in the framework of linear 
cosmological perturbation theory for the SM Higgs inflation models that we consider. A 
review of chaotic inflation formalism in modifled gravitational theories can be found in [40]. 



2.1 Higgs driven ^-inflation 

The general action for these models in the Jordan frame is given by [40, 41]: 
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(2.1) 



where U{(f)) is a general coefficient of the Ricci scalar, R, giving rise to the non-minimal 
coupling, G{(f)) is the general coefficient of kinetic energy and V{(1)) is the general potential. 
The generalized f/(0)R gravity theory in Equation (2.1) includes diverse cases of coupling. 
For generally coupled scalar field C/ = (7 + n^C.cfP'), G{(j)) = 1 and 7 and C are constants. The 
non-minimally coupled scalar field is the case with 7 = 1 while the conformal coupled scalar 
field is the case with 7 = 1 and (^ = 1/6. The conformal transformation for the action given 
in Equation (2.1) can be achieved by defining the Einstein frame line element: 

"^ -^' -dP + a^{i)dx^) , Q = 2k^U{(I)), U{(p) = {I + k^C<P^) , (2.2) 



ds^ = n 



where the quantities in the Einstein frame are marked by caret. From the above equations 
we get: 

\/?7da, di=VUdt, 



da = Vflda , di 

and scalar potential V{(j)) in the Einstein frame is given by: 

1 y(0) 



V{cP) 



4^4 [/2(</,) 



(2.3) 



(2.4) 



The non-minimal coupling to the gravitational field introduces a modification of the Higgs 
field propagator by a suppresion factor s{(j)) defined as [4, 49]: 
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(2.5) 



that made the kinetic energy in the Einstein frame canonical with respect to the the new 

scalar field (p. 

Background equations. The Friedmann equation in the Einstein frame reads as [42, 44]: 



H^ = - 



where: 
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a di 
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(2.6) 

(2.7) 
(2.8) 



Equations (2.7) and (2.8) are enough to compute the background field evolution in the Ein- 
stein frame once the field equations in the Jordan frame are known. 
The Jordan-frame field equations from action (2.1) are obtained as [43]: 
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(2.9) 



that in the slow-roh approximation ( \<p/(j)\ <^ H and |0^| ^ ^('?^) ) can be written as: 
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^ -^7^^772T^('A)> (2-11) 



^ 1 + k2C02(1 + 6C) ^ ^ 



Cosmological perturbations. Introducing the following quantities 

7.^™^^)^, %^V/¥TI74, (2.13) 

the power spectrum of the curvature perturbations in the Einstein frame is given by [43, 44]: 

where r = — 1/[(1 + eaH] is the conformal time. 

The spectral index of the curvature perturbations obtained as usual as ns = l+dlnPs{k)/dlnk 

is: 

ns - 1 = 3 - 2% = 3 - V47S + 1 . (2.15) 

The power spectrum of the tensor perturbations in the Einstein frame is obtained from 
Equation (2.14) making the following replacements: 

<5s -^ <5t = 1 , is -^1T = .-, , -^2 ' Ss-^ St = -777— = . (2.16) 

[i- -r e) 2HQt 

The spectral index of the tensor perturbations n-r and the tensor-to-scalar ratio R are then 
given by: 



riT 



-V^^ ..si(£|))\ ,.n) 



2.2 Higgs driven G-inflation 

In this section we consider that inflation is driven by the SM Higgs potential with a kinetic 
term modified by a function G{(p, X) such that the standard part of the Lagrangian given in 
Equation (1.4) remains much larger than the Galileon interaction term, taking the form [31]: 

£g(</', X) = P{4>, X) - G{^, X)a^ , where \P{X, </.)[ ~ V{4>) » \G{X, (P)a(P\ . (2.18) 



Background equations. The energy density and pressure of the homogeneous and isotropic 
background field described by the FRW metric are [32, 37]: 



p = 2PxX -P + 6G^xH(t>X - G^^X , 
p = P-{G,4, + G,x'^)X, 



(2.19) 
(2.20) 



where H = a/ a is the Hubble parameter. The gravitational field equations are then given 
by: 



and the scalar field equation of motion reads as: 

Px(<A + ^H4>) + 2PxxX4> + 2P^x^X - P 



P, 



2G^X(f>X 



+6G,x[{HX) + 3H^X] - 2G,00X + 6G,xxHXX = . 



(2.21) 



(2.22) 



Since we are interesting in fluctuations generated during inflation, we assume a background 
for which the quantities \(f)/{H(j))\ and \G^(f)(j)/{GH)\ are much smaller than unity. Under 
these assumptions, the first slow-roll parameter can be written as [34]: 



H 



XPx + 3GxH(t>X 



and the sound speed of scalar perturbations is given by: 

2 _ dxp T 



<1, 



where: 



T = Px + ^HG,x, 

g = Px + 2XPxx + &H^{G^x + XG^xx) 



(2.23) 



(2.24) 



(2.25) 
(2.26) 



Cosmological perturbations. The CMB temperature anisotropy directly relates to the statis- 
tical properties of the primordial curvature perturbations TZ on large scales: 



n = -^Q, 



(2.27) 



where (j) is the background value of the field and Q is the perturbation in the unitary gauge, 
<j){x,t) = (p{t), with S(j) = 0. The second order action in this gauge at leading order in the 
slow-roll approximation reads as: 
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(2.28) 



where r is the conformal time. One should note that in the case of G-inflation the uni- 
tary gauge does not coincide with the comoving gauge [31] as in the case of K-inflation for 
which the scalar field Lagrangian is L^ = P{X,(j)) [47]. Consequently, there are two inde- 
pendent solutions of the perturbation equation resulting from the second order action (2.28). 
Introducing the following parameters 



G,x(I^X 



e + u , 



(2.29) 



the power spectra for Q and TZ are obtained as [34]: 

Pq= , 2^2 - ' P'^ = t4i?2—~^ (2-30) 

which are evaluated at the time of sound horizon exit Csk = aH. Defining the additional 
slow-roll parameters: 

the spectral index of the curvature perturbations is obtained as: 

n,-l= '^\' =-2e-fi-s. 2.32 

dhiK 

The power spectrum and the spectral index of tensor perturbations are given as usual in the 
form: 

2H^ 

P. = -^, nr = -2e. (2.33) 

Thus, the tensor-to-scalar ratio R is: 

R=^ = -8csinT-2i^), (2.34) 

that highlights the difference between G-inflation model and K-inflation model for which 

R = —SCsTlT- 

Hereafter we will consider Higgs potential driven G-inflation corresponding to the fol- 
lowing setup [32]: 

G{^, X) = -g{(t>)X , g{(t>) = 4 > (2-35) 

where w is a mass parameter assumed to be positive. Under the slow-roll conditions, the 
energy density is dominated by the potential energy and the standard part of the kinetic term 
in the Lagrangian given in equation (2.18) takes the form \P{(I),X)\ ~ ^(0)- The speed of 
sound corresponding to this setup is c^ ~ 2/3 and the Friedmann equation and the equation 
of motion of the scalar field read as: 

12 ^ ' ^ 3<7(0) ^ ' 

2.3 Higgs driven E-inflation 

In the recently proposed SM Higgs inflation model (E-inflation) with non-minimal derivative 
coupling to gravity [25], the tree-level SM Higgs Lagrangian takes the form: 

Ce{X, ^) = -\ {g^" - w'G'^'') d.cpd.cp - V{<P) , G'^- = i?^- - 1^^- , (2.37) 



where R is the Ricci scalar, G^"^ is the Einstein tensor, w is a mass parameter and V{(f)) is 

the Higgs potential. 

Background equations. The Friedmann equation and the field equation of motion are given 

by: 



6 



l\l + 9H'w') + ^c^' 



dt 



a-'<p{l + 3H^w^) 



3\ j,3 



-a^X(j) 



Asking the solutions to obey the following constraints: 
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< 1 , 9F^u; 



2„,,2 i2 
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^0^ i?» ^ 
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<.m[ 



Equations (2.38) can be written as: 

12 ^ 



SHw'^k'^c 



(2.38) 



(2.39) 



(2.40) 



Cosmological perturbations. The power spectrum of the curvature perturbations evaluated at 
the time of sound horizon exit Cgk = aH is obtained as solution of the perturbation equation 
resulting from the canonically normalized second order action [26]: 



Pn 



K^H^ 
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47r2 12u;2k2</,2 V 2 



where the sound speed of the curvature perturbations is 



0<ci 



3 - ISw^K^d)"^ 



3 + 18u;2k2(/-2 
Introducing the first two slow-roll parameters: 



< 1. 
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K W (p 



Hc^ ^ 2 ' 

to leading order in slow-roll approximation P-ji can then be written as: 



47r2 8e' 

and the spectral index of the curvature perturbations is 

_ dln^T^ _ dlniJ dine 
dln/c dluA; dln/c ' 

and by using Equations (2.43) and (2.44) can be written as: 

ns = l-2e-2T] = l-5e. 
The power spectrum of the tensor perturbations is obtained as: 



Pt 



47r202c3 



(2.41) 

(2.42) 

(2.43) 
(2.44) 

(2.45) 

(2.46) 
(2.47) 

(2.48) 



where and the sound speed of the gravitational waves Cg are: 

fi2 = 1 _ iu;2^2^2 ^ c^ ~ 1 + w'^K^^)'^ . (2.49) 

2 ^ 

To the leading order in uP'k^cJ)^ we have fi ~ 1 and Cg ~ 1. The spectral index of the tensor 
perturbations is 

d InPr 

and the tensor-to-scalar ratio R is given by: 

R=-^ = l2w'^ti^^'^ = Se. (2.51) 

3 Quantum corrections to the Higgs potential 

The quantum corrections due to the interaction effects of the SM particles with Higgs boson 
through quantum loops modify Higgs scalar potential from classical expression in both Jordan 
and Einstein frames, taking the RG forms: 

^^'^ = ~^^*^ ' ^^'^ = T^ (i + .2e(t)02(,))2 ' * = 1- (;^J ' (3-1) 

where the scaling variable t normalizes the Higgs field and all the running couplings to the 
top quark mass scale rriTop = 171.3 GeV [48]. We compute the various i-dependent running 
constants by integrating the RG /3-functions as compiled in the Appendix. For each case, 
the i-dependent running constants are obtained as: 

^(*)=r T^^d*'' Y = {g,g',gs,ytAX], (3.2) 

Jt=o 1 + 7("j 

where 7 is the Higgs field anomalous dimension given by Equation (6.7) from Appendix. 
At t = 0, which corresponds to the top quark mass scale tut op-, the Higgs quadratic coupling 
A(0) and the top Yukawa coupling yi(0) are determined by the corresponding pole masses 
and the vacuum expectation value v = 246.22 GeV as: 



2v' 



A(0) = ttI [1 + 2Ah("^^)] ' ym = - ^[l + AT(mT)], (3.3) 



where Aff(?7T,//) and A j- (7715^) are the corrections to Higgs and top quark mass respectively, 
computed following the scheme from the Appendix of Espinosa et al.(2008). The initial 
values of the gauge coupling constants at t = are given by: 

^ ^^^ - 0.03344 , ^-^ = 0.01027 , ^^ = 0.1071 . (3.4) 



47r 47r 47r 

The couplings g and g' are obtained by integrating RG flow equations from their values at 
M^, while gs is calculated numerically [49]. The value of non- minimally coupling constant C 
at the beginning of inflation is determined such that the calculated value of the amplitude 
of the curvature density perturbations agrees with the measured value [50]. 



4 RESULTS 

4.1 The CMB Angular Power Spectra 

We obtain the CMB temperature anisotropy and polarization power spectra by integrating 
the coupled Friedniann equation and equation of motion of the scalar field with respect to 
the conformal time corresponding to the SM Higgs inflation variants, as presented in the 
previous section. We take wavenumbers in the range 5 x 10^^ — 5 Mpc~^ needed by CAMB 
Boltzmann code [51, 52] to numerically derive the CMB angular power spectra and a Hubble 
radius crossing scale fc* = Csk = 0.002Mpc~^ and impose that the electroweak vacuum 
expectation value v =246.22 GeV is the true minimum of the Higgs potential at any energy 
scale (A(t) >0). 

The value of the Higgs scalar field at the beginning of inflation (f)_\f determines the value 
of the scaling parameter t = \n{(j)j\f/mTop) at this time. 

For the case of ^-inflation the value of (j)j\f relates to the quantum scale of inflation cpj and 
to the duration of inflation expressed in units of e- folding number M through [2] : 

Ai = A (2^4 + (^2 + ^,2)2 _ ig^4) _ g^ ^ (42) 

where the inflationary anomalous scaling parameter Ai involves a special combination of 
quantum corrected coupling constants [6, 49]. 

The value of the field at the end of G-inflation is obtained from Equations (2.23) and (2.35) 
asking the solution to obey the condition e = 1 at the end of G-inflation: 

<^ew = 23/4a-V8^^i/2 . (4.3) 

The number of e-folds till the end of G-inflation is then given by: 

r<i>e.r,d fj I \l/2 1 

from which the value of the field evaluated N e-folds before the end of G-inflation is: 

(l)j^ = (16AA + Sf/^X^^'^Kw^/"^ . (4.5) 



In a similar way, the value of the field at the end of E-inflation is obtained from Equations 
(2.40) and (2.43) by asking the solution to obey the condition e = 1 at the end of E-inflation: 

32 

From the evaluation of the number of e-foldings M during E-infiation: 

f'4>end U 1 

■^ = / T^"^ = H^^^'^' (-^AT - 0U , (4.7) 
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we obtain the value of the field at M e-folds before the end of E-inflation as: 

32 



^'^ = ZJZ^ (SA/" + 1) • (4.8) 



As the inflationary observables are evaluated at the epoch of horizon-crossing quantified 
by the number of e-foldings M before the end of the infiation at which our present Hubble 
scale equalled the Hubble scale during infiation, the uncertainties in the determination of 
M translates into theoretical errors in determination of the infiationary observables [53, 54]. 
Assuming that the ratio of the entropy per comoving interval today to that after reheating 
is negligible, the main uncertainty in the determination of M is given by the uncertainty in 
the determination of the reheating temperature Tr after inflation. 

The dynamics of the Higgs field during G-inflation show that the effect of Galileon-like 
interaction during reheating is very small and can be safely ignored because g = (j)/w^ 
{(j) <C w) term is suppressed around the minimum of the potential [32] . 

In the case of E-inflation the constraint wH ^ 1 suppresses the canonically normalized field 
(f) ~ (wH)^^ that becomes ineffective during reheating. 
Recent studies of the reheating after (^-infiation estimates T^ in the range [5, 55]: 



A 



1/4 



3.4 X 10^3 GeV < T^ < ) 1.1 x 10^^ GeV , (4.9) 

\0.25 / 

that translates into a negligible variation of the number of e-foldings with the Higgs mass 
(AA^ ~ 0.1). The number of e-foldings at Hubble radius crossing scale fc* is related to T^ 
through [38]: 

\ 1/3 / ^3\ 1/3 



AA = ln 



Pend J \9*^r 



HX 



(4.10) 



where g^ ~ 106.5 for the SM, (70 — 2 and Tq is the present value of the photon temperature. 
Taking the variation range of Tr given in Equation (4.9) and A/" = 59 e-foldings at k^,, in the 
view of WMAP7-FSN-FBA0 normahzation at this scale, we obtain 57.8 < TV < 60.2 for the 
number of e-foldings at which A:* exits the horizon. 

Throughout, we will consider M varying in the above range for all SM Higgs inflation variants. 
In this way we conservatively included the same error AM ~ 1.17 to account for possible 
uncertainties in the theoretical estimate of Tj. . 

For each wavenumber k our code integrates the /3-functions of the t-dependent running con- 
stant couplings in the observational infiationary window imposing that k grows monotonically 
to the wavenumber A:*, eliminating at the same time those models violating the condition for 
inflation e = -H/H^ < 1. 

4.2 Markov Chain Monte Carlo (MCMC) Analysis 

We use MCMC technique to reconstruct the Higgs fleld potential and to derive constraints 
on the infiationary observables and the Higgs boson mass from the following datasets. 
The WMAP 7- year data [45, 46] complemented with geometric probes from the Type la 
supernovae (SN) distance-redshift relation and the baryon acoustic oscillations (BAG). The 
SN distance-redshift relation has been studied in detail in the recent unified analysis of the 
published heterogeneous SN data sets - the Union Compilation08 [56, 57]. The BAG in 
the distribution of galaxies are extracted from Two Degree Field Galaxy Redshift Survey 
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(2DFGRS) the Sloan Digital Sky Surveys Data Release 7 [58]. The CMB, SN and BAG 
data (WMAP7+SN+BA0) are combined by multiplying the likelihoods. We use these mea- 
surements especially because we are testing models deviating from the standard Friedmann 
expansion. These datasets properly enables us to account for any shift of the CMB angular 
diameter distance and of the expansion rate of the universe. 

The likelihood probabilities are evaluated by using the public packages COSMOMC and 
CAMB [51, 52] modified to include the formalisms for the SM Higgs driven G-inflation 
and E-inflation as described in the previous sections. 

The fiducial model is the ACDM standard cosmological model described by a set of parame- 
ters receiving uniform priors. For the case of G-inflation and E-inflation we used the following 
set of parameters: 

{Qbh'^ , Oc/l^ , 9s, T , rnuiggs , w , N] , 

where: fJ^/i^ is the physical baryon density, J7c/i^ is the physical dark matter density, 9s is 

the ratio of the sound horizon distance to the angular diameter distance, r is the reionization 

optical depth, ruHiggs is the Higgs boson pole mass, w is the mass parameter and A/" is the 

number of e-foldings at the Hubble radius crossing. 

For the case of ^-inflation we used the same set of input parameters except for the mass 

parameter w. For this case we use as input the amplitude of the curvature perturbations 

Ag ~ A/C^ that fix the value of the non-minimally coupling constant C at A;^, for a given value 

of the Higgs boson mass. 

For each inflation model we run 64 Monte Carlo Markov chains, imposing for each case the 

Gelman &: Rubin convergence criterion [59]. 

Figure 1 and Figure 2 show the dependence of the spectral index of the curvature 
perturbations and of the tensor-to-scalar ratio respectively on the Higgs boson mass, as 
obtained from the fit of the SM Higgs inflation variants to the WMAP7-hSN-hBA0 data set. 
Table 1 presents the mean values and the lower and upper intervals of the input and derived 
parameters as obtained from their posterior distributions corresponding to each model. 

We find that the confidence regions of ns and rriHiggs for G-inflation and E-inflation 
models are correlated as A^g ~ A increases for heavier Higgs boson. For (^-inflation model we 
find the same confidence regions anticorelated. Since for this case we fixed the non-minimal 
coupling constant ^ such that the amplitude of the curvature perturbations ^| ~ A/^^ is at 
the observed value, this anticorrelation reflects the fact that in the inflationary domain the 
running of C dominates over the running of A. Figure 1 clearly shows that at 2-a level the 
variation ranges of ng obtained for SM Higgs inflation variants do not overlap. However, 
to discriminate among different models the measurement of the Higgs boson mass and an 
improved detection of ns are required. 

The connection between the physics of inflation and the low energy physics expected to be 
proved by the collider experiments depends on the choice of UV-completion of each theory. 
For ^-inflation and i?-inflation this introduces an uncertainty in the determination of Higgs 
inflaton mass [14, 60]. For the case of G-inflation it has recently been argued that the choice 
of the Lagrangian in the form L = G{X, (p) leads to second order equations of motion for any 
choice of G{X,(j)) and therefore to unitarity evolution as quantum field theory [61]. 
Therefore, if PLANCK satellite, designed to measure ns to 2-a accuracy of 4 x 10~^, will find 
ns smaller then 0.97, while LHC will find a Higgs mass significantly smaller then 151 GeV, 
the (^-inflation and E'-infiation models would be ruled out. On the other hand, if PLANCK 
will find ns higher then 0.97 while LHC will find the Higgs mass significantly higher then 
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Figure 1. Two-dimensional marginalized joint probability distributions (68% and 95% CL) showing 
the dependence of the spectral index of the curvature perturbations ns on the Higgs inflaton mass 
rriHiggs as obtained from the fit of the SM Higgs inflation variants with totop=171.3 GeV and u=246.22 
GeV to the WMAP7-I-SN+BA0 data set. All parameters are computed at the Hubble crossing scale 
A:*=0.002Mpc-i. 



151 GeV the G-inflation model would be ruled out. One should note that the actual value 
of ns measured by the WMAP team at 68% CL is ns = 0.968 ± 0.012 [45]. 

Figure 2 presents the equivalent dependences for R. The joint confidence regions of R 
and niHiggs are uncorrelated for all SM Higgs inflation variants. These can be attributed to 
a smaller contribution of the tensor modes to the primordial density perturbations when the 
Higgs mass is increased. We also find a striking difference among the predicted intervals for 
R by different models. While in the case of i5^- inflation and (^-inflation the predicted values 
for R are very low, the G-inflation model predicts R values large enough to be detected by 
the PLANCK satellite. Also for this case we conclude that the measurement by LHC of 
the Higgs boson mass and the detection by PLANCK of the tensorial modes are required 
to discriminate among the SM Higgs inflation variants. If PLANCK will detect a value of 
i? ~ 0.1 while LHC will find a Higgs mass significantly smaller then 151 GeV, the (^-inflation 
and E'-inflation models would be ruled out. 
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Figure 2. Two-dimensional marginalized joint probability distributions (68% and 95% CL) show- 
ing the dependence of the tensor-to-scalar ratio R on the Higgs inflaton mass muiggs as ob- 
tained from the fit of the SM Higgs inflation variants with ?titop=171.3 GeV and w=246.22 GeV 
to the WMAP7+SN+BA0 data set. All parameters arc computed at the Hubble crossing scale 
fc*:=0.002Mpc-^ 



Figure 3 presents the dependence of the tensor-to-scalar ratio on the spectral index of 
tensor perturbations (the consistency relation) as obtained from the fit of the SM Higgs infla- 
tion variants to the WMAP7-|-SN-|-BA0 data set. Figure 4 clearly shows that the consistency 
relations are unique enough to distinguish the SM Higgs inflation variants. 

Finally, we address the question of the violation of perturbative unitarity for E-inflation 
model. Figure 4 presents the dependence of the scale at which unitarity is violated in this 
model, A ~ {2H'^k)^''^, on the value of the background field (j) and on the comoving curvature 
perturbations TZ ~ 12H^, as obtained from the MCMC analisys. We find that A is smaller 
than both \/R and (p during E-inflation, confirming that this model suffers from unitarity 
problems [27], contrary to the original claims [26]. 



5 CONCLUSIONS 

Although the SM Higgs inflaton non-minimally coupled to gravity can account for inflation 
(("-inflation), other variants of Higgs inflation within the SM based on scalar field theories 
with non-canonical kinetic term such as Galileon-like kinetic term (G-inflation) and kinetic 
term with non-minimal derivative coupling to the Einstein tensor (E-inflation) have been pro- 
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Figure 3. Two-dimensional marginalized joint probability distributions (68 % and 95% CL) showing 
the degeneracy between the spectral index of tensor perturbations nx and the tensor-to-scalar ratio R 
obtained from the fit of SM Higgs inflation models with 77iTop=171.3 GeV and ■y=246.22 GeV to the 
WMAP7-I-SN-I-BA0 dataset. For each case the dashed line shows the theoretical consistency relation. 
All parameters arc computed at the Hubble crossing scale A:H,=0.002Mpc~^. 

posed. In this paper we have analyzed the possibility to distinguish ^-inflation, G-inflation 
and E-inflation models through their predictions for ns and R as a function of mniggs as 
well as through their predictions for the consistency relations. 

In order to ensure consistent results, we have studied the SM Higgs inflation variants consider- 
ing the full dynamics of the background and perturbations of the Higgs field during inflation 
and compute the RG improved Higgs potential by including the two-loop /3-functions for 
the SU(2) X SU(1) gauge couplings, the SU(3) strong coupling, the top Yukawa coupling, 
the Higgs quadratic coupling and the Higgs field anomalous dimension. Assuming that all 
SM Higgs inflation variants are consistent theories, we use the MCMC technique to de- 
rive constraints on the inflationary parameters and the Higgs boson mass from their fit to 
WMAP7-FSN-FBA0 data set. We show that the SM Higgs inflation variants lead to signifi- 
cant constraints on Ug, R and mniggs- 

From the analysis of the confidence regions in ns-rriHiggs plane we conclude that in order to 
discriminate the SM Higgs inflation variants by using these parameters the measurement by 
LHC of ruHiggs and an improved detection of ns are required. If PLANCK satellite, designed 
to measure ns to an accuracy of 4 x 10~^, will find ns smaller then 0.97, while LHC will 
find a Higgs mass significantly smaller then 151 GeV, the (^-inflation and £^-inflation models 
would be ruled out. 
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Figure 4. Two-dimensional marginalized joint probability distributions (68% and 95% CL) showing 
the dependence of the scale at which unitarity is violated A on the value of the background field 
4> (left) and on the comoving curvature perturbations TZ (right) as obtained from the fit of the E- 
inflation model with mTop=171.3 GcV and ■i;=246.22 GeV to the WMAP7+SN+BA0 data set. AU 
parameters are computed at the Hubble crossing scale A:*=0.002Mpc~^. 



We also find a striking difference for the predicted values for R by different models. While 
for i?-inflation and (^-inflation the predicted values for R are low, the G-inflation model pre- 
dicts R values enough large to be detected by the PLANCK satellite. We conclude that 
if PLANCK will detect a value of i? ~ 0.1 while LHC will find a Higgs mass significantly 
smaller then 151 GeV, the ^-inflation and i?-inflation models would be ruled out. 
We also show that the consistency relations are distinct enough to differentiate the SM Higgs 
inflation variants. 
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Table 1. Mean values and the lower and upper intervals (at 68% CL) of parameters obtained 
from the fit of the SM Higgs inflation variants with mTop=171.3 GeV and ■y=246.22 GeV to the 
WMAP7+SN+BA0 data set. All parameters are computed at the Hubble radius crossing /c,=0.002 
Mpc"i. 
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6 Appendix 

In this appendix we collect the SM renormalization group /3-functions at renormalization 
energy scale t = \ii[(l)/mTop) beyond the top quark mass rnxop- 

All /3-functions and the Higgs field anomalous dimension 7 include the Higgs field suppression 
factor s{t) that for the case of ^-inflation is given by [4, 38]: 

For the case of G-inflation and E-inflation we take s{t) = 1. 

The two-loop /3-functions for gauge couplings gi = {g',g,gs} are [7]: 



/5g. = Kgfbi + K^gf 



3 



Y,B,jg]-sit)diy^t 
j=i 



(6.2) 



where k = I/IGtt and 



/ 199/18 9/2 44/3 \ 
6= ((40 + s(t))/6,-(20-s(t))/6, -7), B=i 3/2 35/6 12 , 

V 11/6 9/2 -26/ 

d* = (17/6,3/2,2). (6.3) 
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For the top Yukawa coupling yt, the two-loop /3-function is given by [4]: 

9 



Pyt = Kyt 



-\9^ - ^9" - 85' + ls{t)yl 



+ i^^yt 



23 4 3 2/2, 1187 /4 „ 2 2 
-^9 -^9 9 + -7j:^9 + 95- Qs 



+ ^49' 9l - 108<7! + f ^5^ + ^5'^ + 365.^ ) s{t)yl + 6 i-2s\t)yt - 2s\t)ylX + sHt)X'' 



9 



V 16 



16 



(6.4) 



The two- loop /3-function for the Higgs coupling A is [4, 38]: 



Px = K 



+ K^ 



/2\2 



I8s\t) + 6)X' - Gyf + ^ (2g' + {g' + g'^Y) + {-9g' - 3g'' + Uyf) X 

^ (915/ - 289/ff'2 _ 559^2^4 _ 379^/6) + 3(3^(^)^6 _ ^4 /8|;! ^ 33^2 ^ 3^^^^^ 



+ A (-|/ + ^5'<7'' + ^s{t)g'' + 108s\t)g'X + 36s\t)g"X - 312s'{t)X' 
+ y? (-^/ + Y^'a" - '{9" + A (f 5^ + f <7'2 + 8O5? - 144.^(t)A 



(6.5) 



The /3-function for non-minimal coupling ^ is given by [38]: 



I3( = ^(C + i ] {Gil + s{t))X + 6y^ - y^ 



;9 



(6.6) 



Finally, the two- loop Higgs field anomalous dimension 7 is given by [4]: 

"271 , 



7 = —K 



+ K^ 



4 4 

/45 



M 



85 



32 



. 9 2/2 431 ,4 

■^ - 16^ ^ - -W'^'^' 



27 



I ¥^ + ^9" + 20<7,^ ) yf + -s{t)yt - Qs\t)X 



(6.7) 
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